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I address the question of what can be learned from the observation of the diffuse supernova 
neutrino flux in the precision phase, at next generation detectors of Megaton scale. An analytical 
study of the spectrum of the diffuse flux shows that, above realistic detection thresholds of 10 MeV 
or higher, the spectrum essentially reflects the exponential-times-polynomial structure of the original 
neutrino spectrum at the emission point. There is only a weak (tens of per cent) dependence on the 
power P describing the growth of the supernova rate with the redshift. Different original neutrino 
spectra correspond to large differences in the observed spectrum of events at a water Cerenkov 
detector: for typical supernova rates, the ratio of the numbers of events in the first and second 
energy bins (of 5 MeV width) varies in the interval 1.5 - 4.3 for pure water (energy threshold 
18 MeV) and in the range 1 - 2.5 for water with Gadolinium (10 MeV threshold). In the first 
case discrimination would be difficult due to the large errors associated with background. With 
Gadolinium, instead, the reduction of the total error down to 10 — 20% level would allow spectral 
sensitivity, with a dramatic improvement of precision with respect to the SN1987A data. Even in 
this latter case, for typical neutrino luminosity the dependence on 13 is below sensitivity, so that it 
can be safely neglected in data analysis. 



I. INTRODUCTION 

Almost twenty years after the detection of neutrinos 
from SN1987A, there is a hunger for new neutrino data 
from supernovae, needed to advance our understanding 
of the physics of core collapse and precious to improve 
our knowledge of neutrino physics. 

Considering how rare supernovae are in our immediate 
galactic neighborhood [H, it is likely that the next 
supernova neutrino data will come to us not in the form 
of a high statistics signal from a nearby star, but as a 
diffuse flux of neutrinos from all the supernovae in the 
sky, each of them too far to give a significant signal. 

There are reasons to believe that the detection of this 
diffuse supernova neutrino flux (DSNi/F) may soon be a 
reality. Indeed, the SuperKamiokande (SK) experiment 
has already obtained an upper bound on this flux within 
an order of magnitude of the theoretical predictions: 1.2 
electron antineutrinos {ve) s~"'^cm~^ Q (see also (1, 
for other, less constraining, experimental results). The 
sensitivity of SK will increase substantially in the next 
5-10 years with the Gadolinium addition, GADZOOKS 
0. In this new configuration, the rate of events due to 
the DSNz/F could reach 0{1) events/year, which implies 
the possibility to obtain the very first detection of the 
diffuse flux. If the supernova rate is precisely known, 
GADZOOKS could also give information on the neutrino 
spectrum at least for part of the parameter space, with 
sensitivity comparable to that of the SN1987A data 

After SK and GADZOOKS, a new generation of detec- 
tors with larger volumes will become operational. Some 
projects will require ~ 1 Megaton water volumes (Hy- 
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perKamiokande, UNO, MEMPHYS H, [M [H), while 
other ideas involve a 50 kiloton mass of liquid scintilla- 
tor (LENA^fllJlll) and 100 kiloton-scale Hquid Argon 
detectors [l4 Il5j|. For values of the DSNz^F that are 
considered typical, these new detectors are expected to 
accumulate a statistics of several tens or even hundreds of 
events from the DSNz/F in the space of few years running 
time. 

It is clear, thus, that the advent of these new detectors 
will make it possible to go beyond the phase of discovery 
and start a phase of detailed study of the diffuse flux. The 
number of events, their energy spectrum, and their spa- 
tial and temporal distributions will be analyzed in depth 
to study the neutrino emission from a supernova, the lo- 
cal and cosmological supernova rate, neutrino oscillations 
and possible exotic particle physics affecting neutrinos. 

The question of what will be possible to learn from 
the DSNz/F data in the precision phase, and how well, 
has received only limited attention. A initial study by 
Yuksel, Ando and Beacom Q considered specific exam- 
ples with focus on the low statistics of GADZOOKS. The 
main result was that GADZOOKS has a chance to test 
the neutrino spectrum if the flux is close to the current 
SuperKamiokande limit. The status of the art on what 
data from the DSNi^F will tell us is well summarized in 
the work by Fogli et al. [l^, where a detailed study 
of the background is presented for detection with pure 
water and with water with Gadolinium. 

In this paper I expand on previous works by discussing 
in detail what can be learned from the energy spectrum 
of the data from the diffuse flux, with focus on Mega- 
ton scale detectors. The energy distribution of events 
depends on the spectrum of the neutrinos emitted by an 
individual supernova and therefore its study would be 
a crucial test of models of neutrino transport inside the 
collapsed matter of a supernova, as well as of neutrino os- 
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dilations inside the star. The observed energy spectrum 
also depends on the distribution of supernovae with the 
redshift, i.e. on the dependence of the supernova rate 
with the distance (or, equivalently, time from the Big 
Bang). Indeed, the larger the fraction of supernovae at 
cosmological distance, the softer the DSN:/F spectrum 
is, due to the effect of redshift of energy. To get infor- 
mation on the evolution of the supernova rate with the 
cosmological time from the diffuse neutrinos would be of 
extreme interest. 

Here I present an analytical study of the DSNj/F spec- 
trum as well as numerical results for the spectrum of v^- 
induced events at a Megaton water Cerenkov detector. 
The first offers insight on the latter. Of both the DSNi/F 
spectrum and the observed signal, I examine the depen- 
dencies on the power of evolution of the supernova rate 
with the distance and on the average energy and spec- 
tral shape of the neutrino flux at the source. I also give 
typical numbers of events and discuss the sensitivity of 
specific observables to the neutrino spectrum, considering 
the errors associated with the background and possible 
uncertainties on the supernova rate normalization. 

The paper opens with generalities (Sec. |TT|, followed 
by the analytical study of the DSNi/F spectrum in sec. 
mil The observable signal is discussed in Sec. |TVl Sum- 
mary and discussion close the paper in Sec. [V] 



II. GENERALITIES 

A core collapse supernova is an extremely powerful 
neutrino source, releasing about 3 • 10^^ ergs of energy 
within ^10 seconds in neutrinos and antineutrinos of 
all flavors. The spectrum of these neutrinos is approx- 
imately thermal, with average energies in the range of 
10 — 20 MeV and with the muon and tau species being 
harder than the electron ones due to their weaker (neu- 
tral current only) coupling to matter. Inside the star, the 
neutrinos undergo either partial or total flavor conversion 
depending on the mixing angle ^13 and on the mass hier- 
archy^ ordering) of the neutrino mass spectrum (see e.g., 
[I3i[3)- For definiteness, here I focus on the Ve species, 
as it is the most relevant for detection (see Sec. IIV[) . 

The energy spectrum of v^, as they exit the star can be 
described as |19| : 

dN_ ^ (l + a)^+°L /^y^-(i+„)B/Bo ^^ 
dE r(l + a)£;o' \E^) ' W 

where Eq ^9 — 20 MeV is the average energy and L 
the (time-integrated) luminosity in Pg, L ^ 5 ■ 10^^ ergs, 
a is a parameter describing the shape of the spectrum, 
a ~ 2 — 5 [19], with larger a corresponding to narrower 
spectrum. The form ^ is justified by a number of argu- 
ments. First, it is a good fit to numerical calculations of 
neutrino transport before flavor conversion. Given this, 
Eq. ID) also describes the Dg spectrum after conversion if 
the conversion is total, which happens for inverted mass 



hierarchy and sin^ ^13 ^ 10~^ [13, EBIj or if the Pe and 
v^tVt original spectra are very similar. This is favored 
by recent finding that the De and P^, Vr average energies 
may differ only by 10 - 20% fl9]. If the conversion of 
Vf, is partial and the hierarchy between the electron and 
non-electron original spectra is strong, Eq. ([1]) is still an 
accurate description of the neutrino spectrum above its 
average energy, where the flux is dominated by the harder 
component from oscillated muon and tau neutrinos p^ . 
This condition is met above the realistic detection thresh- 
olds of ~ 10 — 20 MeV energy, since the DSN:/F average 
energy is typically smaller than ~ 5 — 6 MeV (see fig. [T]). 

I have estimated the error associated with approximat- 
ing a composite spectrum with the form ([T|), in the case 
of maximal permutation of the neutrino fluxes, which is 
realized (for antineutrinos, of interest here) for the nor- 
mal mass hierarch y, o r for the inverted hierarchy with 
sin^ ■ 10~^ nil- I assumed a difference in the 

Ve and Vpj^jUr average energies up to 30% and up to a 
factor of two difference in the luminosities of the electron 
and non electron flavors. The result is that the error is 
smaller than 6% above 10 MeV of energy, if Eq. ([T]) is 
used with L and Eq being the luminosity and the aver- 
age energy of the composite spectrum, and a chosen to 
reproduce the width of the composite spectrum. 

Finally, the form ^ is supported, as phenomenological 
description, by the fact that it fits well the data from 
SN1987A [33. 

Given the neutrino output of an individual supernova, 
one has to sum over the supernova population of the uni- 
verse. This is described by the cosmological rate of su- 
pernovae, Rsn{z)^ defined as the number of supernovae 
in the unit time per unit of comoving volume at redshift 
z. Its value today (z = 0) is estimated in the range of 
i?sw(0) 10^"* Mpc^'^yr^^, with uncertainty of a factor 
of two or so [13, Hlj due to finite statistics and uncertain 
systematic errors due to dust extinction. Theory '22'| and 
data indicate (see e.g. [2l|) that after the beginning of 
star formation the supernova rate (SNR from here on) 
was first somewhat constant in time, and then decreased 
after the redshift dropped below z ^ 1. It can be de- 
scribed phcnomcnologically with a broken power law: 

Rsn{z) = RsNmi + zf forz<l 

= RsNm^ for z > 1 , (2) 

where j3 is in the range /3 ~ 2 — 5 from the statistical 
analysis of supernova observational data only [20| . It 
becomes more constrained around /3 ~ 3 (best fit value 
(3 = 3.28) if the more precise (but only indirectly related 
to the SNR, see [20l . [23j ) measurements of the star for- 
mation rate are used instead [2l|. Here I will consider 
the wider interval /? 2 — 5 to be conservative. 

Using the supernova rate Rsn and the individual neu- 
trino flux, dN/dE, one finds the diffuse flux, differential 
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m energy: 



HE) - jT 



Rsn{z 



dN{E') 



dz 



dE' 



v/an(l + 2)3 + f^A 

(see e.g. [23j). Here E and E' = {\+z)E are the neutrino 
energy at Earth and at the production point; ~ 0.3 
and f2A ~ 0.7 are the fraction of the cosmic energy den- 
sity in matter and dark energy respectively; c is the speed 
of Hght and TJo is the Hubble constant. The maximum 
redshift Zmax accounts for the expected decrease of the 
SNR with z beyond z ~ 5 or so. The dependence of the 
DSNi/F above realistic detection thresholds on Zmax is 
weak due to the effect of redshift. I will use Zmax — 5 
from here on. 

Throughout the paper, the expression ([3|) will be used 
for the DSNi/F. This formula is valid in the continuum 
limit, and only if individual variations between different 
stars are negligible. This last statement has not been 
checked in detail; it is however supported by the results 
of numerical simulations 'ss']. 

From Eqs. H]), (l2|), and ([3|) one sees that the spectrum 
of the DSNi^F depends on three variables: a, Eq and 
/3, i.e. on the spectrum of the neutrinos emitted by an 
individual star ("original neutrino spectrum" from here 
on) and on the power of growth of the SNR with z. Fig. 
[1] shows the DSNt/F spectrum from Eq. ([3]) for different 
original neutrino spectra; in the figure L, i?5Ar(0) and 
[3 are held fixed while a and E^ are varied (see caption 
for details) to show how the slope of the spectrum can 
change. Specifically, the spectrum is steeper for smaller 
i?o and larger a, as it will be discussed in detail in Sec. 

nn 



III. THE ENERGY SPECTRUM OF THE 
DIFFUSE FLUX 



What is the energy distribution of the diffuse super- 
nova neutrinos? How is it related to the original neu- 
trino spectrum and to the cosmological distribution of 
supernovae? In the light of what discussed in Sec. [ill 
the answer is only partially intuitive. Due to the red- 
shift of energy, the farther a supernova is from us the 
more its contribution is shifted to lower energy. This 
means that the DSNz/F does not retain the same spec- 
trum that the neutrinos had at emission; but have in- 
stead a distorted - and overall softer - distribution with 
respect to it, which depends on the SNR. One expects the 
largest distortion in the lowest energy part of the spec- 
trum ~ 5 — 10 MeV, where the contribution of high 
redshift supernovae is larger. In contrast, in the high en- 
ergy end of the DSNi^F spectrum, ^ 15 - 20 MeV, the 
flux is largely dominated by stars with negligibly small 
redshift, and therefore it should be closer to the original 
spectrum. 

Let us see this in a more quantitative way, by examin- 
ing the spectrum that results from Eq. ([3]). There, the 
integration can not be done exactly. Still, however, one 
can find useful analytical approximations that help to 
understand the energy dependence of the DSNi^F. Con- 
sidering that current and planned searches of the diffuse 
flux have thresholds of ~ 10 MeV or higher, it is sensible 
to focus on the high end of the spectrum, i.e. on the 
small redshift approximation, z <C 1. In this limit one 
can make the following simplifications: 



neutrinos/Me V/s/cm' 



^=12 MeV, a = 5 
^=15 MeV, a = 2.6 
^=20 MeV, a = 2 




1. The substitution: 



[n,^{i + zf + nK 



-1/2 



(1 + ^) 



(4) 



The two expressions coincide at first order in z un- 
der the condition that -I- f^A = 1, which is fa- 
vored by current cosmological data. 



2. Neglect the contribution of the supernovae beyond 
the z ~ 1 break of the SNR, specifically: 

Rsn{z) = i?5Ar(0)(l z)^ for Z Zmax , 

= for z > z,„ax , (5) 



FIG. 1: Examples of spectra of the diffuse supernova neutrino 
flux for different original neutrino spectra and fixed super- 
nova rate. For the latter, I used the parameters /3 = 3.28, 
Rsn{Q) = 10"'' Mpc"^yr"\ The luminosity L = 5 • 10^^ ergs 



(see Eq. ([T}) was adopted. 



With these approximations, the integral in Eq. ([3]) 
jives the following result: 
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Rsn{0) 



L 



Ho r(2 + a)e2 



a—r}—l 



(r [77 + 1, S/e] - r [ry + 1, (1 + z„,a.)i;/e]) , 
I 



(6) 



where I define ry = a+f3—3ilm/2 and e = i?o/(l+a) for 
brevity. Considering a ^ 2 (Sec. |n|, we liave e ^ Eo/3. 

It is possible to further simplify the result by drop- 
ping the contribution of the upper limit of integration 
(1 + Zmax)E/ e, which at high energy is exponentially 
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suppressed with respect to the term due to the lower 
limit. Furthermore, if 77 is integer, one can also use the 
function in terms of an exponential 
H, and get: 



expression of the T 
times a polynomial 



$e ^ Rsn{0) 



or, explicitly: 



RsNiO) 



L 



E 



c L _E X ^ 



Ha r(2 + a)e2 



+ 77(77-1) 



E 



a-3 



r 



(77 -fc)! 



(7) 



+ 77(7/-l)(77-2) 



(8) 



Eq. ([8]) turns out to be numerically accurate also for 
non-integer 77, provided that the sum is truncated at the 
integer nearest to 77. It also has the advantage of simplic- 
ity, so it can be used to understand the physics. 

Let us pause for a moment and study the result ([8]). 
This expression tells us that the DSNi^F spectrum is sim- 
ilar to the original spectrum in the general structure of an 
exponential times a polynomial part. More specifically, 
the term e^^/^/e^ in Eq. ([5]) is the same as in the original 
spectrum, while the polynomial part is more complicated 
and contains the information on the cosmology (i.e., on 
(3 and on fim), through the quantity 77. The latter enters 
the result only in the coefRcients of the various powers of 
energy, a very weak dependence with respect to the ex- 
ponential decay of the spectrum in E / e. Notice that, in 
agreement with intuition, the dependence on 77 becomes 
weaker for higher energy: indeed, in the limit E ^ e, 
the polynomial part is dominated by the 77-independent 
term {E/e)°'~^. In general, Eq. ([5]) suggest that the ex- 
ponential captures the most of the energy dependence, 
with only slower variations due to the polynomial part. 
Thus, one can consider fitting the DSNi^F spectrum at 
high energy with a simple exponential form, as already 
suggested in \2^: 



$e(0)e- 



(9) 



I find that, expectedly, (E) is numerically close to e, with 
differences smaller than ~ 30%. 

Fig. [5] compares the exact spectrum (see caption for 



details), calculated numerically, with the three approxi- 
mations: the sophisticated one, Eq. the cruder one, 
Eq. (O, and the phenomenological one, For the lat- 
ter, the parameters (E) and $e(0) have been adjusted 
to fit the numerical curve in the interval _E ~ 20 — 30 
MeV. The figure shows that the sophisticated result, Eq. 

underestimates the flux at low energy as a result of 
neglecting z > 1 part of the SNR. Instead, the the cruder 
description, Eq. overestimates the fiux at low energy 
due to dropping the (negative) contribution of the upper 
limit of integration in ([S]), which is equivalent to taking 
the SNR as growing with constant power all the way to 
z — > 00. The simple fitted exponential form lies between 
the two approximations and works surprisingly well down 
to i? '--^ 5 — 6 McV, probably because by fitting (E) one 
effectively reproduces at least part of the contributions 
that are not included in the forms (O and (O, especially 
the z > 1 part of the SNR. 

I have studied the precision of Eqs. ^ and ([7]) in 
detail. First, I fixed /3 = 3.28, and varied the neutrino 
spectral parameters in the intervals a = 2 — 5 and Eq = 
9-18 MeV. I obtained that the fuU resuh, Eqs. ®, is 
precise to better than 12% (40%) a.t E ^ 19.3 MeV {E = 
11.3 MeV). Instead, the cruder expression ([7]) exceeds the 
exact result by up to 40% at E = 19.3 and up to a factor 
of three at = 11.3 MeV. 

The analytical results lose accuracy for larger (3, as ex- 
pected, since larger /3 means larger contribution of super- 
novae at z > 1. With /3 = 5, the sophisticated result ^ 
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FIG. 2: An illustration of the approximate analytical descrip- 
tions of the DSNj/F spectrum in comparison with the exact 
spectrum calculated numerically (details in the legend). The 
lower panel is a close-up of the upper one on linear scale. I 
used the same SNR parameters and luminosity as in fig. [T]and 
the spectral parameters Eo — 15 MeV, a = 2.6. For the fit- 
ted exponential form, Eq. ®, the value of the characteristic 
energy is (E) = 4.57 MeV. 



errs by up to 20% (120%) a.t E = 19.3 MeV (E = 11.3 
MeV). Eq. ([7]) can deviate from the exact result by a 
factor of several. 

The exponential form, Eq. ([9]), with (E) fitted to the 
reproduce the exact result at E' ~ 20 — 30 MeV, is accu- 
rate within 5% everywhere above 20 MeV and deviates 
from the exact spectrum by at most 25% at E = 11.3 
MeV. 

In closing this section, I would like to emphasize that, 
while useful for understanding, the analytical approxi- 
mations discussed here lose validity below 15-20 MeV of 
energy. Thus, to model observables at these low energies, 
it is necessary to calculate the neutrino flux numerically 
using Eq. ([T]). I do so in the reminder of this work. 



IV. DETECTING THE DIFFUSE FLUX: 
OBSERVABLES 

A. Inverse beta decay: signal and background 

I now turn to discussing the spectrum of events induced 
by the diffuse flux in a detector. In particular, I consider 
events due to inverse beta decay, +p n + e~^, since 
they largely dominate the signal in water and liquid scin- 
tillator. 

The detected positrons from inverse beta decay carry 
information on the incoming spectrum. This is thanks 
to the fact that, except for small effects due to nucleon 
recoil, the neutrino and positron energy differ simply 
by the mass difference of neutron and proton: E^+ = 
E — 1.29 MeV. Thus, the positron energy spectrum is 
roughly given by the neutrino spectrum distorted by the 
energy dependence of the detection cross section, which 
goes approximately as E^. 

To calculate the spectrum of the observed positrons, 
I model the diffuse i?e flux according to Eq. ([T]), with 
Zmax = 5. In all cases, I will use the typical values 
Rsn{0) = 10"* Mpc-^r-i and L = 5 ■ 10^^ ergs, and 
vary a, Eq and (3 in the intervals that I indicated as typ- 
ical in Sec. [TTl I use the detection cross section from ref. 

For definiteness, results are given for a water Cerenkov 
detector of 0.45 Mt fiducial volume (20 times the volume 
of SK) without and with Gadolinium, after 4 years of 
operation. I consider two different energy thresholds, 18 
MeV and 10 MeV in positron energy, characteristic of 
water only and water with Gadolinium (water-|-Gd here- 
after) respectively. I adopt the same energy resolution 
and efficiency (~93% in the energy range relevant here 
[1^) as SK [2^, |23|. The set ups considered here are the 
same as those in ref. [l6j (exce pt f or a minor difference 
in the volume: 0.4 Mt in ref. [3 and 0.45 Mt in this 
work), so that a direct comparison is possible. 

For a meaningful interpretation of the results, it is nec- 
essary to take into account the background carefully. As 
described in a number of publications (e.g., [23| and ref- 
erences therein), there are four sources of background 
that in water result truly ineliminable. These are: spal- 
lation neutrons, reactor neutrinos, atmospheric neutrinos 
and invisible muons. The first one is dealt with by set- 
ting the energy threshold at 18 MeV of positron energy. 
Above this cut the spallation and reactor backgrounds 
are practically zero. The remaining backgrounds, atmo- 
spheric neutrinos and invisible muons, can only be in- 
cluded in the statistical analysis, as done in ^3,], so that 
the signal and these backgrounds can be separated only 
on statistical basis. Of these, the invisible muons back- 
ground is larger than the DSNj/F signal by a factor of 
several [si, fl6|. 

The addition of Gadolinium to the water will allow 
to discriminate inverse beta decay events from events of 
other type, thanks to the detection in coincidence of the 
positron and of the neutron capture on Gd [7]. With 
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Gd the energy threshold can be lowered to 10 MeV of 
positron energy, motivated by the presence of the reac- 
tor i/g's below it. It is estimated that in the water+Gd 
configuration the invisible muon background will be re- 
duced by a factor of '^5 with respect to pure water 
This reduction will bring it down to be smaller or com- 
parable to the signal above the 10 MeV threshold. 

Here I model the background following ref. [l^. As 
done there, the events due to background will not be 
shown, but will be included in the calculation of the sta- 
tistical error. More explicitly, I consider a scenario in 
which from the total of the data one subtracts the back- 
ground using a Monte Carlo prediction for it. The error 
on the subtracted signal is then given by the statistical 
error on the total (signal plus background) number of 
events: a = yjNsig + N^ckg- 

The results of this work can be easily rescaled to repro- 
duce the inverse beta decay signal at other detectors (up 
to differences due to different energy resolution). The 
case of the proposed 50 kt liquid scintillator detector 
LENA [l^, [l^ is particularly interesting. Since scintil- 
lator distinguishes inverse beta decay through the de- 
tection of the positron and neutron in coincidence, its 
performance will be similar to that of a water-|-Gd de- 
tector, up to a rescaling to account for the different vol- 
ume. I estimate that, in terms of number of events, 1 (4) 
year(s) of running time of a 0.45 Mt water -|-Gd detec- 
tor with the efficiency given above corresponds to ^7.4 
(^30) years running time of LENA with 100% efficiency. 
Equivalently, for the same running time, LENA will have 
a number of events smaller by a factor 0.13. 



B. The positron energy spectrum 

Figures [3] and [4] give the spectra of events in 5 MeV 
bins for the water only configuration and the water+Gd 
one. For visual convenience, in each figure all the his- 
tograms have been normalized at the same number of 
events, N=60 for water only and N=150 for watcr+Gd. 
These are the numbers obtained with a = 2.6, Eq — 15 
MeV and (3 = 3.28. Each figure has two panels: in the 
upper one I fix the cosmology and vary the neutrino spec- 
tral parameters, a and Eq^ while in the lower panel the 
spectral parameters are fixed and (3 is varied (see fig- 
ure captions for details). The error bars (omitted in the 
lower panels for simplicity) represent the 1 a total (sig- 
nal+background) statistical errors. 

Let us first comment on the water only case. Fig. [3l 
It appears immediately that the spectrum changes sub- 
stantially when the spectral parameters are varied. The 
behavior agrees with what expected from Sec. IIIII the 
spectrum falls more rapidly for the smaller e, according 
to the exponential dependence in the DSNj/F spectrum. 
So, a first conclusion is that in principle the signal in- 
duced by the diffuse flux is very sensitive to the original 
neutrino spectrum. 

In spite of this, however, different DSNt/F spectra - 



and therefore different original spectra - will be difficult 
to distinguish one from the other with high statistical 
confidence. This appears from the error bars plotted in 
fig- El One can see that the largest variation of the num- 
bers of events in the first bin is within 2 a statistical er- 
ror, and the variations are smaller in the other energy 
bins. As already anticipated in [iGj . this means that 
with typical statistics a Megaton water detector would 
be able to give indication or evidence of the DSNj/F, but 
not to give precise spectral information about it. With 
higher DSNz/F flux - within what is allowed by the cur- 
rent SK limit - the spectral sensitivity could be mod- 
erately better. For example, with twice as many signal 
events, the variation of signal between the two extreme 
cases in the first energy bin would be of about ~ 3cr. No- 
tice that the total statistical error is dominated by the 
invisible muon background and therefore would increase 
only weakly with the increase of the signal. 

The lower panel of Fig. [3] confirms the expectation of 
little dependence on the SNR parameter j3. Even with 
the rather extreme variations of (3 used in the figure, 
the number of events varies by at most ~ 15% in the 
bins where the statistics is higher. Considering the large 
errors, I conclude that very little or no information on (3 
could be extracted from DSN:/F data. 

With water-l-Gd the potential to study the DSNi/F 
spectrum is dramatically better: indeed, the statistical 
errors are smaller thanks to the improved background 
reduction, and the signal is enhanced by the lower en- 
ergy threshold. Fig. 2] illustrates this. It appears that 
the two extreme cases differ by more than 3(t in the first 
bin. This means that the possibility of spectrum dis- 
crimination is concrete, even though the differences in 
the other bins are not as significant. As in the case of 
pure water, the dependence of the observed spectrum on 
[3 is weak. The difference between the number of events 
in the first bin for the two extreme scenarios is at the 
level of ~ Icr, increasing to ^ (2 — 2.5)(t for doubled or 
tripled signal. Thus, to discriminate between different 
values of [3 would be difhcult, even in the ideal case of 
perfectly known original spectrum. 

The points made above can be illustrated more quan- 
titatively using the ratio r of the number of events in the 
first and second energy bin. Taking bins of 5 MeV width 
this ratio is: 

_ jV(18<ge+/MeV<23) 
^ 7V(23 < E^+ /MeV < 28) ' ^ ' 

for the water only case, and 

_ N{IQ < E^+ /MeV < 15) 

~ 7V(15 < £:e+/MeV < 20) ' ' 

for water-|-Gd, in terms of the positron energy E,,-!- . To 
limit oneself to the numbers of events in the first two bins 
is overall more convenient with respect to including the 
data in the higher energy bins. The latter add little to 
the signal but add substantially to the background, and 
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FIG. 3: The spectra of the observed positrons from inverse 
beta decay for the pure water setup, for fixed SNR (upper 
panel, /3 — 3.28 used) and for fixed original neutrino spectrum 
(lower panel, Eq = 15 MeV, a = 2.6 used). In all cases the 
numbers of events have been normalized to A'^ = 60. 



therefore they end up washing out the effect of the signal 
by increasing the error. 

For pure water and fixed SNR (Fig. [3l upper panel) I 
find that r varies in the range r ~ 1.5 — 4.3. With the 
statistics used in Fig. [31 the error on r is larger than 
^ 100%; this confirms the conclusion that the sensitivity 
to the neutrino spectrum is limited. For fixed original 
spectrum and varying /3, the variation of r is minimal, 
much smaller than the error. 

For the water+Gd configuration in the case of fixed 
SNR (Fig. m upper panel) I have r ~ 1 — 2.5. The error 
is ~ 20 — 30%, meaning that at least the extreme cases 



FIG. 4: Same as fig. |3]for water+Gd. The numbers of events 
have been normalized to A'' = 150 in all cases. 



should be distinguishable for typical luminosity of the 
DSNi^F. For fixed original spectrum and varying f3 I get 
r ~ 1.2 — 1.6, which is smaller than the error unless the 
DSNi/F is close to the SK upper limit. 



C. Number of events 

To complement the study of the DSN^^F spectrum - 
which is the main focus of this paper - here I give a 
brief discussion of the event rates and their information 
content. 

As it is clear from the previous sections, the number 
of events in a detector depends on the original neutrino 
spectrum, on the SNR power /3 and on the overall factors 
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Rsn{0) and L, i.e. the normalization of the SNR and the 
neutrino luminosity in P,, (after oscillations). Once the 
detection of the DSNi^F moves from the discovery to the 
precision phase, the event rate will constrain the product 
^SAr(O) X L and will complement the spectral analysis to 
extract Eq, a and f3. Realistically, these two aspects can 
not be separated: the data will be analyzed by a global fit 
with both the spectral parameters and the overall factors 
as fit variables. Still, for illustration it is useful to discuss 
certain limiting cases where the sensitivity to the overall 
factors and that to the spectral parameters can be treated 
separately. 

If the DSNz/F spectrum is reconstructed very precisely 
by the spectral analysis, the event rate will serve mainly 
to constrain Rsn{0) x L. This can happen for the 
water+Gd setup, where the statistical errors are small 
enough to allow spectral reconstruction. High precision 
would be obtained if the DSNz/F signal is particularly lu- 
minous, typically two or more times larger than in Fig. 
m In the ideal case of perfectly known spectrum, the er- 
ror on the Rsn{0) x L inferred from the data would be 
just the statistical error, of the order of tens of per cent. 
This would be an enormous improvement with respect 
to the current information from supernova observations 
(which constrain Rsn(O)) and from the neutrino data 
from SN1987A (see e.g. 

It should be considered that by the time the DSNt/F 
is detected, the SNR could be better known from new 
supernova surveys like JWST [H] and SNAP since 
these are expected to see thousands of core collapse su- 
pernovae while undertaking their primary mission of ob- 
serving type la supernovae. From these surveys we are 
likely to obtain a precise measurement of (3, while some 
systematic uncertainty on Rsn{0) will be left due to ex- 
tinction effects. In this situation, the measured /? will 
be used as input in the DSNi^F spectral analysis, thus 
increasing the precision of reconstruction of Eq and a. 
The astrophysical measurement of Rsn (0) and the infor- 
mation on Rsn{0) X L from the neutrino event rate will 
be combined to constrain L and reduce the systematic 
error on Rsn{0)- If the systematics due to extinction 
is well understood and taken into account, the DSNi/F 
event rate will essentially measure L. 

It is possible that the product Rsn{0) x L and the 
power P will be much better known than the original 
neutrino spectrum. This scenario could be realized if 
RsNiO) and f3 are precisely measured by supernova sur- 
veys (plus precise treatment of systematic effects) and L 
is determined by improved numerical simulations of core 
collapse, perhaps combined with astrophysical data on 
supernovae such as energy of the explosion, mass of the 
remnant, etc.. 

In this case the event rate will serve to constrain the 
spectral parameters in combination with the energy dis- 
tribution of the data. In fig. [5] I briefly illustrate this 
scenario in the ideal situation in which L, Rsn{0) and 
/3 are known with negligible error (see caption for de- 
tails). The figure shows isocontours of the ratio r and 



of the event rate N in the plane a — Eq, both for water 
only and for water with Gadolinium. It appears that the 
event rate varies by a factor of several depending on the 
spectral parameters: from ~ 10 to ~ 140 (~ 70 — 230) 
events in 4 years for water only (water plus Gadolinium) 
[3^ . The figure also illustrates a possible measurement 
of N and r, with central values and Icr statistical errors 
on these two quantities [37] . If only r is considered - for 
robustness against uncertainties on Rsn{0) x L - then 
for water only a major portion of the parameter space is 
allowed within Icr, and the entire space is within 3a. In- 
stead, for water -|-Gd the la allowed region is much more 
restricted, and the 3a region excludes the edges of the pa- 
rameter space; for example, the measurement of r = 1.2 
would exclude the point Eq = 12 MeV and a = 4.5. In- 
cluding the information on the number of events reduces 
the allowed region substantially: while for pure water the 
whole a ~ Eq plane is allowed within 3a, for water -|-Gd 
the potential for exclusion is better. For the specific case 
in the figure, a measurement of = 200 would exclude 
Eq = 12 MeV regardless of the value of a. 

Notice that the region allowed by A^ is contained inside 
the one allowed by r, meaning that the number of events 
alone is sufficient to constrain the original neutrino spec- 
trum, if L, Rsn{0) and /3 are known very precisely. In 
presence of errors on these quantities, which is proba- 
bly more realistic, the ratio r will be more constraining 
instead. 



V. SUMMARY AND DISCUSSION 

It is a concrete possibility that the next supernova neu- 
trino data will be from the diffuse flux. If so, these data 
will be of interest to improve on our knowledge of the 
supernova neutrino spectrum with respect to what is al- 
ready known from SN1987A. Furthermore, the DSNz^F 
data will be analyzed to test the cosmological rate of su- 
pernovae. 

Here I have addressed the question of the sensitivity 
of the DSNz/F to these quantities, both at the level of 
the i>e flux at earth and of the spectrum of inverse beta 
decay events in a detector. 

One flrst answer is that the spectrum of the diffuse 
De at Earth above realistic detection thresholds (10 MeV 
or higher) essentially reflects the neutrino spectrum at 
the source (after oscillations), with only a weak (tens of 
per cent) dependence on the power f3 of evolution of the 
supernova rate with the redshift. Indeed, with respect 
to the original spectrum, the DSNz^F spectrum retains 
the same structure of an exponential times a polynomial 
part. The exponential term is identical to the one in the 
original spectrum, and dominates the energy dependence 
of the DSNz/F spectrum. The information on the cosmol- 
ogy is in the coefficients of the polynomial, and therefore 
the cosmological parameters influence the spectrum more 
weakly. These features, that were studied in detail here, 
are in agreement with the general intuition about the 
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FIG. 5: Solid (red) lines: isocontours of the number of events 
in 4 years running time, A^, in the space a — Eq. Long dashed 
(black) lines: isocontours of the ratio r of the numbers events 
in the first and second bin (Eqs. (|1U|I and ([TT}). In each 
panel I show a possible measurement of A'' and r, with central 
values (diamond) and la statistical errors (regions within the 
short dashed lines, marked with arrows; the wider region is 
the error on r). The SNR parameters and luminosity are as 
in fig. m 

contribution of more distant supernovae being weak at 
high energy due to redshift. 

At the level of experiment, a better back- 
ground subtraction with respect to pure water 
(e.g, by the addition of Gadolinium to the water 
of UNO/HyperKamiokande/MEMPHYS or by using a 
large scintillator detector (LENA)) would be needed to 
move from the discovery phase to the precision phase. 

What major advances will come from the precision 
phase depends on the situation in the field at that time. 
A likely scenario is that there will be new, precise mea- 
surements of the supernova rate. These will result in 



a precise determination of /3, while leaving some uncer- 
tainty on the overall scale of the SNR, Rsn{0), due to 
effects of extinction. In this case, the DSN^F data will 
mainly improve our knowledge of the original neutrino 
spectrum, thus adding precious information to what is 
already known from SN1987A. For typical luminosity of 
the signal, after 4 years of operation a 0.45 Mt detector 
in the water-|-Gd configuration would have hundreds of 
signal events, and total (signal-|-background) statistical 
error of 10 — 30% in the first two energy bins. This 
should make it possible to discriminate at least between 
the extreme cases of neutrino spectrum. 

In addition to the spectrum, the observation of the 
DSNi^F will allow to measure the product of the neutrino 
luminosity and of Rsn(0). This measurement will be 
combined with direct supernova observations to reduce 
the uncertainty on Rsn{0). 

In the exciting case in which data from a new galactic 
supernova become available before the DSNi/F is mea- 
sured, the study of the DSNj/F will take a completely 
different direction. The high statistics data from the 
galactic supernova will give the neutrino spectrum and 
luminosity with very good precision. These and the value 
of P from supernova surveys will be used as input in the 
analysis of the DSN:/F to extract information on the less 
precise parameter i?5jv(0). 

Depending on the precision on the DSNi^F signal, a 
number of interesting tests could be performed. One 
of them is to look for differences between the value of 
-RsAf(O) estimated in astronomy and the one obtained 
from neutrinos: the comparison could reveal a number of 
"failed" supernovae: objects that undergo core collapse 
with the emission of neutrinos without a visible explo- 
sion. It would also be interesting to look for differences 
between the neutrino spectrum inferred from an individ- 
ual galactic supernova and the one extracted from the 
DSNi^F: such differences will be a measure of the vari- 
ation of the neutrino output from one supernova to an- 
other. Differences could be expected in consideration of 
the fact that the DSNj/F receives the largest contribu- 
tion from low mass supernovae, M ~ 8 — IOA/q (simply 
because the initial mass function goes roughly like the 
power —2 of the star's mass, see e.g. [2l| and references 
therein) , whose neutrino output could differ from that of 
a supernova with a larger mass progenitor. 

In summary, the present study shows the relevance of 
studying the energy spectrum of the diffuse supernova 
neutrino flux, mainly as a tool to reconstruct the 
original neutrino spectrum and therefore to test the 
physics of core collapse that influences the neutrino 
spectrum formation inside the star. It also motivates 
the experimental effort to reduce the background in the 
relevant energy window in order to achieve the necessary 
precision for spectral sensitivity. 

I acknowledge support from the INT-SCiDAC grant 
number DE-FC02-01ER41187. 
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